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We obtain a generalized Christoffel-Darboux (GCD) formula for skew-orthogonal polynomials. 
Using this, we present an alternative derivation of the level density and two-point function for 
Gaussian orthogonal ensembles and Gaussian symplectic ensembles of random matrices. 
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\^ ' Random matrices have found applications in different branches of physics mainly due to the 'universality' in their 
, correlation function under certain scaling limits. In this context, although unitary ensembles of random matrices, 
corresponding to systems with broken time-reversal symmetry (for example, a mesoscopic conductor in the presence 
[ of a magnetic field) have been extensively studied much less is known about Orthogonal Ensembles (OE) and 

■ Symplectic Ensembles (SE) of random matrices. 

^ [ In his phenomenal paper Q Dyson stressed the importance of a good understanding of the skcw-orthogonal pply- 
'"^ i nomials (SOP) to study the OE and SE of random matrices. Some progress has been made in this rega rd [H, 13, M to 
OO ' develop the theory of SOP. Alternative approaches were taken by various authors [gIF^ 1^ HolllTl lia. Il3 | to study 
. these ensembles. For example, Deift and Gioev 0, have recently used the Widom's representation |l6j | to prove 
Universality for a wide class of polynomial potentials of the OE and SE of random matrices. In all these work, the 
^ , authors have used the well known properties of orthogonal polynomials to study these ensembles. Here we use the 
' more elegant SOP, evolving naturally from OE and SE of random matrices. 

In this paper we obtain recursion relations for SOP. Using this we derive the GCD formula. We use them to obtain 
I the level density and two-point correlation function for Gaussian orthogonal ensembles and Gaussian symplectic 
' ensembles of random matrices. 

, We consider ensembles of 2N dimensional matrices H with probability distribution 

O ■ 1 

■ Pp,N{H)dH ^ ^exp[-[TrV{H)]]dH, (1) 



X 



where the parameter (3—1 and 4 corresponds to H real symmetric or quaternion real self dual. (Note that we have 
considered the OE of even dimension. The odd dimension can be easily generalized.) dH is the standard Haar 
measure. 2pN is the so called 'partition function', and is proportional to the product of skew-normalization constants 



/np-l 
exp[-[Try(i/)]]dff - iV! n 5f \ (2) 

3=0 

Here np is called the 'Fermi level' by analogy with a system of fermions. In our case, it takes the value np = 2N 
for P = 1 and 4 respectively. It is assumed that most of the physical quantities (like level density and correlation 

functions) are related to properties of gif ^ around the 'vicinity' 18] of the Fermi level. In this paper, we give rigorous 
justification for such a claim. 

To study different correlations among the eigenvalues of random matrices, we need to study certain kernel functions 
0,1^. For example, the two-point correlation function for /? = 1 and 4 can be expressed in terms of the 2x2 matrix 



i^^y) - I 71) ' Jp) I' (3) 



while the level-density 



p^^\x)^si^^ix,x). (4) 
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Here, the kernel functions are defined as 

2N 

= ^'^^\y)\{^'-^\x), (5) 
D'^ij^{x,y) = m{x)\{^^P\y), (6) 

N 

lf^{x,y) = -^^^\x)\{^^^\y)+5,/-^^^ (7) 

2N 

where 

$(« = (<I>f ...$(f)...)*, $(« = -$(«*z, (8) 
(similarly for \I>('')) are semi-infinite vectors. They are formed by quasi-polynomials 



where 



and 



(9) 



(x) cxp[—V{x)], (10) 



^\?\-) = E^"'"^-'' (11) 

fc=0 



^ = ( ° J )+...+ (12) 



is the SOP of order n. 

. ( 

is a semi- infinite anti-symmetric block-diagonal matrix with = —1 and 

r 

5 is the kronecker delta. The matrix 

[] = diag(l,...,l,0,...,0) (13) 

have 2A^ entries. Finally, we define 

M/l^) (a:) = <i>f ) (x), (a;) ^ j $(1) (2;),(:, - y)dy, (14) 

which satisfy skew-orthonormality relation P, 0, : 

($(«,ci>03))^ / $(f )ci>L^)dx = (15) 



The contom of integration T' is the real axis. Here, one must mention that under suitable assumptions on V{x) that 
ensures convergence of the integral for a given V{x), these SOP are unique up to the addition of a lower even 
order polynomial to the odd ones. 
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For the unitary ensemble, study of correlation function involve similar kernel function, which is calculated '3, Q 
using the well known Christoffel Darboux formula To study the kernel functions arising in OE and SE (|5i6l7f) . 
we derive a GCD formula. 

To this effect, we expand x^'^^'^^x), (<i>^'^)(x)) and {x^^^\x)) in terms of ^^^\x) (and hence introduce the semi- 
infinite matrices Q''^\ -P^^-* and R'^^^ respectively): 



*(4)(a;) ^ ^^(4)^(4) 



(16) 
(17) 
(18) 



where (|18ll is obtained by multiplying the above expansion by e(y — x) and integrating by parts. They satisfy the 
following commutation relations: 

[g(«,p(«] = i, [i?(«,p(^)] = p(«. (19) 

Using {'lJJn\x),^pm\x)) and (xipn^x) , 2lJm\x)) for /9 = 4, and replacing ip^^^{x) by (/i(-^^(x) for /? = 1, we get 

pw ^ _piP)^^ pm ^ (20) 



where dual of a matrix A is defined as 



= -ZA*Z. 



(21) 



However starting with {x(j>i^\x),^pm\x)) and using H19|l . we get 

Q(/3) ^ Q(/5)^ + (F(/3))"\ (22) 

We will use the matrices P^^'' and i?'''^ to obtain the GCD formula. For /3 = 4, using ((SJ, (^Tjl and if^ . we get 

$w*(x) n (y) + *^'^*(-) n ^ n "^^'^ 

2N 2N 2N 2N 

$(4)*(x) n ^ n + 'i'^'^ V)^* n ^ n '^'^'^ (y^ 

2N 2N 

-$(4)*(a;)ZZ[];z[|P$(*)(y) 

2Af 2N 

^^'^^^\x)zzp'zz n ^ n 



2N 2N 



$(*'(.t) 



2N 2N 



$W(y), 



(23) 



where we have used ■^n2iv ^Tl2N ~ ~ Tl2N- Similarly, using lO, ({TUl . ifTTIl and ipn|l we get 

yS^^j^{x,y) - a;S'^^(y,x) 

= $(4)*(a;)P*J]zn'f*'^(y) + $''^*(a^)n^n^'^*'^(2/) ' 



2N 2N 



2N 2N 



¥^\x) 



(24) 



2iVj 
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Finally using ^ and the GCD for the SE (/? = 4) is given by 
Following a similar procedure, the GCD for the OE (/? = 1) is given by 



'5'27v(2^i y) 



y-x 



Here, one might recall, that for orthogonal polynomials, the Christoffel-Darboux sum takes the form 

:(2) 



S'^{x,y) = <^\x)\[^iy)^ ^'^^^^^^^-^^^y^ 



N 



x-y 



(25) 



(26) 



(27) 



where Q is a tri-diagonal Jacobi matrix, and are the normalized orthogonal quasi-polynomials. Any matrix of the 
form [A, Hjv] only off-diagonal blocks whose size depends on the number of bands above and below the diagonal 
of A. Thus the 'size' of P'^^^ and r!'^^ decides the number of terms around the 'Fermi-level' that will ultimately 
contribute to the correlation. For example, for the orthogonal polynomials, the Christoffel-Darboux has only two 
terms. 

Having derived the GCD for arbitrary weight, we will study ensembles with 



d+l 

V{x) = ^ 



1=1 



Ul I 

-^' 



(28) 



where V{x) is a polynomial of order d+l. is called the deformation parameter 0|. We will show that for such 
ensembles, the matrix P^^^ and R'^^^^ are finite band matrices, i.e. have finite number of bands below and above the 
principal diagonal. 



For /3 = 4, 



while for (3—1, 



^ ^ ^n.rn't'rri (•^) 

[-v'ix)4^H^) + ^l'U^) + - 

-E(^'(Q)W™ (^)+'^n-l(^) + •■ 



^ [(l^nHy) - y)dy 

[-E(^'W))n,™'^™^(y)] <^-y)dy 

+VWi(x) + .... 



(29) 



This gives 



Similarly, for (3 — A, we have 



— lower. 



(30) 



(31) 



= hxy'(a;)</>W(x) + n44)(^) + ... 

= -(EQ"™[(^'(^?^'^))""^i''(^)) + 



(32) 
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while for /3 = 1, we get 



^ \y(t>n Ky) - y)dy 
-yV (y)(l)^n\y) ^{x-y)dy 

+ {n + l)tljll\x) + . . . 

m.l 



(33) 



Thus we get 



i?(/3)+Q(/3)y'(g(/5)) 



= lower- 



(34) 



where 'lower' denotes a strictly lower triangular matrix and 'lower'|_ a lower triangular matrix with the principal 
diagonal. Since Q'^^'^ has only one band above the diagonal, Eas. H31l34|l confirm that P*^''-' and i?*^^-' has d and c? + 1 
bands above and below (since they are anti-self dual) the principal diagonal. However, unlike /3 = 2, Q^''^ for (3—1 
and 4 is not a finite band matrix. 



Q'^^^ can be calculated in terms of the normalization constant gif'^ and the coefficients of the polynomials c 
For SOP, with g^^j^ = 52w+i' have 



„(fe./9) 



9j+i 



S-1 



(13) ' ^3,J {j,l3) 



„(i+i,/3) 



Finally using Eas. (P3|) . and the asymptotic results for these SOP 0, we present an alternative derivation of the 
level-density and the 'two-point' function for the Gaussian orthogonal and symplectic ensembles. For the Gaussian 
Ensembles, with d—1, (jSHJ and if^ give 



and 



^2n{x^ y) 



ra;p(4) _r(4) 1^(4) , r(4) ,(4) _n(4) ,(4) 
1-^^(1,2) ^(1,2)JV'(0,2) + -"'(0,2)'>"(1,2) ^(1,3) ''^(0,3) 

X — y 



(35) 



respectively, where 



5'2^(x, y) — 



ly^(l,2) -"■(l,2)J''^(0,2) -"-(0,2)'^(l,2) -"-(1,3) ''^(0,3) 



a(4) 



y~x 



(4) 



(and similarly for V-'lj fc) ) and 



4(/3) ^ .(/3) 

{j,k) — ^2N+jaN+k- 



For /? = 4, using gf^ = {2n + l)\n^'^2'^'' and c^,"'^^ = (^2)"^" ' 5], we have = 1/2^ and q[^J2) 

From Ea. (pTT|l and with U2 = 2 and mi = 0, we get 



P, 



(4) 



-2gj^;2) = ~AN/V2, 

;?(4) _ „9n(4) r)(4) - _n 

-"•(0,2) — ^'^(0, 1)^(1, 2) ~ 



(4) 
(1,3) 



^^(1,2)^(2,3) ~ ■ 



(36) 
(37) 

(38) 
V2A^. 

(39) 
(40) 
(41) 
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For large TV, we use the asymptotic results (with Gaussian weight) for the SOP Q: 



where 



4%) 



sin [f^'^Hn,e)] 



ni/^VTT sin 6' 



1 



{An) 



1/4 



(42) 
(43) 



(n + 3/4) (sin 26*- 261) + ^, 
2 / /9(n, x)dT + 

-\/2;r 4 



(44) 



and a; = (2r7, + 3/2)"'^^^ cos 0. For a given x, 9 depends on n; for example On — On^^i — ±(2ntan6'„) ^. Then with 
y = X + Ax and 9 = 9„, where 0„ = (02ri^2n+i) , and expanding in 9 and N, we get from the first term in Ea. H35|l : 



^p(4) r ,(4) , 



1 COS 61 . f ,^df^^HN,9)\ . „^ 
sin A6> ' sin 26*. 



961 



(45) 



The second and third terms in Ea. (|35|l give 

'-(0,2) ^''^(1,2) 



n(4) r,(4) i _ _ p(4) r,(4) 
-^(0,2)1'>"(1,2)J ^ ^(1,3) ^''^(0,3) 

COS 26* 



47r sin 6* 



■ sin A0 



,df^ '\N,9) 
89 



(46) 



We know that the odd SOP is arbitrary to the addition of a lower even order polynomial. This choice cancels the 



term R 



^(i'2)\-'^\fi2)^- Collecting all the terms, we get 

sin[2^(2iV-a;2)Ax] 



^2Ni^^ y) 



2ttAx 

Taking Ax 0, we get the famous 'semi-circle', while 



|x| < V27V. 



S'2'^(x, y) sin[27rAxS'^'' (x, x)] sin27r7- 



■'2N\ 

:(4) 



JL 

(4), 



27rr 



(47) 



(48) 



where r — AxS'2^(x,x), gives the universal sine- kernel in the bulk of the spectrum. 



For ^ = 1, we have from 0] g^^^ = (27i)!7ri/222" and c 
'N. For M2 = 1 and ui = 0, Eqs.(|SU and (P^l gives 



(2n,l) 



„(2n+l,l) 
-2n+l 



(1) 
(1,2) 



= 2^", which gives Qjo'i) ~ 



P, 



(1) 

(1.2) 



(1) 

(1.2) 



^(na) ^ ^(0,1)^(1,2) ~ 
n(i) _ n(i) 

^(1,3) ^ ^(1,2)^(2,3) 



= -TV. 



(49) 
(50) 
(51) 



For large iV, we use the asymptotic results for the SOP 0: 

.;i) _ 



1 

'2ni/4 



cos 



[/«(", 0)] 



(52) 
(53) 
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where 



/(i)(n,e') = (n + l/4)(sin26'-26l) + ^, 



= / p{n,x)dx + —, (54) 
J-VIn 4 

and a; = (4n + 1)^^^ cos0. Writing y — x + Ax, and expanding in and in we get 



,(1) _ sin[v/(47V-a;2)Aa:] 



Taking Ax — > 0, we get the level density, while 



(56) 



S2^^{x,x) '!tAxS2^^{x,x) 

where r = AxS^^{x, x), gives the universal sine-kernel in the bulk of the spectrum. 

For general d the correlation function corresponding to a weight with single support can be obtained using the 
asymptotic results for the SOP 18] . However, one needs to understand in greater detail the structure of the finite- 
band matrices and hence the matrix Q^'^-' to come up with a proof. 

In conclusion, the unitary ensembles of random matrices, which involve the orthogonal polynomials have been well 
studied in recent years. In contrast, barring a few specific weights, nothing much is known about the OE and SE. 
This is mainly due to the hurdles created by the SOP. 

In this paper, we have made some progress in understanding some of the basic properties of these SOP. In this 
context, we would like to emphasize that the GCD formula, derived in this paper, can at best be considered as the 
first step for a systematic study of the OE and SE of random matrices. One still needs to develop the theory further 
to come to an equal footing with the unitary ensemble of random matrices. For example, one would like to understand 
in greater details the asymptotic behavior of these SOP to study different correlations for a larger family of OE 
and SE. We would also like to point out the similarity in the GCD formula for (3=1 and 4 with the interchange of $ 
with ^E*. This may be useful in proving the duality between these two ensembles. We believe that these SOP satisfy 
a. d X d differential system, which can be used to formulate a Riemann-Hilbert problem for these matrix models. We 
wish to come back to a few of these questions in a later publication. 

I am grateful to Bertrand Eynard for useful discussions. 
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